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We construct the strong-coupling series in 4d simplicial quantum gravity up to volume 38. It is used to calculate 
estimates for the string susceptibility exponent 7 for various modifications of the theory. It provides a very efficient 
way to get a first view of the phase structure of the models. 
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1. INTRODUCTION 

Euclidean quantum gravity is formally defined 
by the partition function 



Z = 



Vol(diff) 

where the Einstein- Hilbert action 
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Discretizing the theory on an ensemble of four- 
dimensional manifolds T, consisting of N4 4- 
simplexes glued together along the faces, and as- 
suming a metric where all link lengths are equal 
to one, leads to simplicial quantum gravity: 
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. K 2 N 2 (T)-K i N i (T) 



(3) 



Here K4 corresponds to the cosmological constant 
A, K2 to the inverse Newton's constant G, and N2 
is the number of 2-simplcxes (triangles). C(T) 
is the symmetry factor of triangulation T, the 
number of equivalent labelings of the vertexes. 
Eq. ([|) is related to the canonical (fixed volume) 
partition function through 
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Although it is not proven analytically, one as- 
sumes that Zjv 4 is exponentially bounded in A/4. 
Moreover, in our analysis we assume that 



7(k 2 



+ * + (5) 



when N4 — > 00, where 7 determines the critical 
behavior of the grand-canonical partition func- 
tion: Zqc ~ (/i c — K4) 2 " 7 as K4 — > /i c . 

Simulating this model one observes two phases; 
a branched polymer phase for K2>1.3, and a 
crumpled phase for K2<1.3. In the branched 
polymer phase 7 = o> whereas in the crumpled 
phase the sub-leading correction is exponential, 
rather than Eq. (^); this corresponding formally 
to 7 = — 00. Separating the two phases is a first 
order phase transition, hence no interesting criti- 
cal behavior is observed. For that one may need 
to modify the simple action in Eq. (JJj) above. 
However, as investigating every modification re- 
quires extensive numerical simulations, it is im- 
portant to have some analytical guidance. One 
such is the strong-coupling expansion. 

We have used this expansion to investigate two 
modifications of the model, described in Ref. [Q, 
namely the interaction of gravity with gauge mat- 
ter fields, and a modified measure. 



2. STRONG-COUPLING EXPANSION 

The strong-coupling expansion consist of ex- 
plicit construction of the smallest configurations 
and their symmetry factors, and a calculation of 
the additional weights corresponding to a given 
model. In two dimensions matrix models allow a 
recursive construction of the series; in four dimen- 
sions, though, such methods are not available, but 
an alternative procedure based on numerical sim- 
ulations has been developed [0: 
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(a) Using Monte Carlo simulations we identify all 
distinct triangulations, up to a given volume. To 
identify the triangulations, we introduce a hash 
function /(T); this function has to be compli- 
cated enough to distinguish between combinatori- 
ally different triangulations, while simple enough 
for a repeated calculation in the MC simulations. 

Only triangulations that can not be reduced by 
applying one of the volume decreasing geometric 
moves have to be considered in the MC. All others 
are systematically constructed from smaller ones. 
This reduces the MC efforts substantially as, in 
practice, 99.9% of the triangulations turn out to 
be reducible. 

(b) For each triangulation we calculate the corre- 
sponding symmetry factor C (T) by comparing all 
possible permutations of the vertex labels (only 
permutations of vertexes of same order have to be 
considered). Comparing the C(T)'s to the rela- 
tive frequency with which different triangulations 
appear in the MC simulations provides a check on 
the correct identification of the triangulations. 

Using this procedure we have identified all four- 
dimensional triangulations up to volume N4, = 38, 
all in all f. 477. 713 distinct triangulations. This 
gives the first 16 terms in the strong-coupling 
series, however, due to different finite-size cor- 
rections these terms split into two distinct series 
which have to be analyzed separately. 

The strong-coupling series is analyzed using an 
appropriate series extrapolation method, e.g. the 
ratio method. Assuming the asymptotic behavior 
Eq. (||) this yields the critical coupling fj, c and the 
exponent 7. 

3. RESULTS 

Using the strong-coupling expansion it is easy 
to explore a large set of modified models of 4d 
simplicial gravity. Here we will give results just 
for two such models: 
(a) Modified measure: 
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where o(tj) is the order of the j'th triangle, i.e. 
to how many 4-simplexes it is attached. 
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Figure 1. Strong-coupling estimates of 7 for both 
a modified measure (solid line) and / copies of 
vector fields (dashed line). The circles are values 
of 7 from MC with modified measure. 



(b) / copies of vector fields: 



Zv(T) = fl[dA(l) 



3 -S[A(i)] 



(7) 



S(A(l)) = 5>(Uc)^(U) + A(l bc ) + A(l ca )} 2 

tabc 

Here A(l) are non-compact U(l) gauge fields re- 
siding on the links and the sum is over all trian- 
gles. The prime indicates that the zero modes of 
the gauge field are not integrated over. 

The weights corresponding to both models are 
calculated explicitly for every configuration T. 
The series so obtained are analyzed with the ra- 
tio method, assuming the asymptotic behavior 
Eq. (|^). In Fig. 1 we shown an example of 7, for 
K 2 > «|> f° r both the models as we vary (3 and 
/ respectively. For consistency we have verified 
that successive approximations in the analysis of 
the series converge, what is plotted is the result 
of the highest approximation. These results are 
in excellent agreement with values of 7 measured 
in MC simulations of Eq. @ for N 4 = 4000. Fur- 
thermore, the two models yield the same variation 
of 7 provided one takes 



P 2 4 



(8) 
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Figure 2. Variations in 7 versus k 2 for / copies Figure 4. The first corrections term in Eq. (5). 

of vector fields coupled to 4d-gravity. 



i 

1 

\ 
\ 


Branched 

v y=l/2 d H =2 


Crumpled 

Y = — 00 

d H=°° 








Crinkled 

y<0 d H >2 



Figure 3. A schematic phase diagram of simplical 
gravity coupled to either a measure term, or to 
vector fields. 



This universality appears in a variety of "reason- 
able" modifications of the model Eq. (^), all of 
which give the same phase structure (modulo a 
trivial re-scaling of parameters) Q . 

Fig. 2 shows 7 as we vary k 2 i for different num- 
ber of vector fields, given by the strong-coupling 
series. The values of 7 become unstable for small 
K2; this is consistent with the fact that in the cor- 
responding MC simulations we observe a phase 
transition to a crumpled phase. This transition 
is evident in the series extrapolation as succes- 
sive approximations fail to converge both in the 
crumpled phase and in the critical region. 



From these and further investigations we con- 
clude that phase diagrams of the two models are 
quite similar, as depicted schematically in Fig. 3. 
In addition to the crumpled and branched poly- 
mer phases, a new crinkled phase appears, for suf- 
ficiently strong coupling to either matter or the 
measure term, characterized by a finite 7 < 0. 

As K2 — > 00, the finite-size effects become, 
rather surprisingly, less important in the pure 
gravity model. This is related to the dominance 
of a particular class of manifolds, stacked spheres, 
in this limit ||. However, as a coupling to a mea- 
sure term is added the finite-size effects increase 
again, as shown in Fig. 4 were we plot the the first 
correction term to Eq. (|^) for K2 — 00 and varying 
(3. The increased finite-size correction coincided 
with the appearance of the crinkled phase. 

The examples given above show that the 
strong-coupling series gives a good qualitative 
view of the phase structure of modified models of 
4d simplicial gravity. It provides us with a power- 
ful tool and should be used before starting exten- 
sive computer simulations of a particular model. 
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